All photons are equal but some photons are more equal than others 
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Photons indistinguishability is studied from a pragmatic, rather than speculative, point of view. 
As consequence of a detailed analysis of the coincidence probability distribution in a Hong-Ou- 
Mandel experiment, we introduce the concept of rate of distinguishability of photons. For photons 
prepared in a given, but arbitrary, state, this parameter allows one to make a quantitative analysis 
of their indistinguishability. Among other things, our study reveals that coupling between different 
degrees of freedom of the electromagnetic field, plays a crucial role in photons indistinguishability. 
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Introduction. — Although particles indistinguishability 
has been a tenet of quantum mechanics long since, it is 
still a not well understood concept Actually, there is 
a large amount of literature on this subject, mostly laying 
in the realm of "foundations of quantum mechanics" , fre- 
quently nearby the edge between physics and philosophy 
[2|-l5j]. However, besides being a central conceptual issue 
of modern physics, particles indistinguishability is also 
the key to experimental tests of fundamental laws of na- 
ture, as the spin-statistic theorem which is a pillar of the 
standard model of particles and forces @, 0| • In quantum 
electrodynamics, particles (photons) indistinguishability 
manifest itself when two equally prepared photons inter- 
fere at the two input ports of a 50/50 beam splitter (BS). 
The joint probability of detection (coincidence probabil- 
ity) at the two outputs is exactly zero. Vice versa, when 
the two photons are prepared in a different manner, the 
coincidence probability raises up to 50%. This effect was 
first demonstrated by Hong et al. 8] and rapidly be- 
came central in a broad range of experiments in quantum 
physics . Very recently, photon indistinguishability 

proved to be of special interest also for quantum infor- 
mation processing [l3, 14 1. Despite such a lively research 
activity, at the present time some misconceptions still re- 
main regarding indistinguishable particles. For example, 
it is a common idea that "[...] quantum mechanical par- 
ticles are too small for one to attach physical labels; they 
don't have enough degrees of freedom to enable one to 
mark each particle differently" fljj . 

In this Letter we demonstrate that, contrary to popular 
belief, particles (specifically: photons) do have enough 
degrees of freedom to permit one to assign different dis- 
tinguishing marks to each particle, as many as the num- 
ber of degrees of freedom (DOFs) of the photons them- 
selves. To this end, we introduce in an operational man- 
ner the rate of distinguishability of photons which ap- 
pears to be a very useful "tool" . Our study shows that 
coupling between different DOFs critically affects pho- 
tons indistinguishability. These results suggest the need 
for replacement of the strong concept of "particles in- 



distinguishability" tout court, with the weaker concept 
of "particles indistinguishability with respect to a given 
degree of freedom" . 

Two-photon interference. — To begin with, let us con- 
sider the two-photon interference experiment sketched in 
Fig. 1. 8]. The two photons impinging upon the 50/50 




FIG. 1. (color online) Two-photon interference at a 50/50 
beam splitter. Two independently prepared photons enters 
the two input ports A and B of the BS. They are eventually 
detected by two distinct detectors placed behind the output 
ports C and D of the BS. The plane of the figure is the plane 
of incidence 



BS are prepared in the product state \ip AB ) = \^ A )\tp B ), 
where \<>p A ) = tf[ip A ]\0) (\^p B ) = tf[ip B ]\0)) denotes the 
single-photon state in arm A (B), with a^[i/'" 4 ] (6^[V'' B ]) 
being the operator that creates one photon in the wave- 
packet mode ip A (ip B ) flrl The normalized func- 
tion ip A - B completely fixes the spectral, polarization 
and spatiotemporal characteristics of the photon enter- 
ing port A, B. Annihilation and creation operators as- 
sociated to orthogonal wave-packet modes, do commute: 
[d[tp], b'lcf)]] — (tp,(f>)5ab, where (ip,^) denotes the scalar 
product in the complex linear space Jzf 9 ij),<f> UM- The 
probability of detecting the two photons at the two out- 
put ports C and D is given by [IS) ]: 



P M [^> B ] = (l-|(^lO| 2 )/2 



(1) 



2 



where {ip A \ip B ) = (ip A ,ip B ). 

Now, assume that the two single-photon wave packet 
modes ip A and ip B are prepared "almost" in the same 
manner, namely ip A = ip an d ip B = (ip + Sif)) /\\ip+Sip\\ 
i. ' + Sip. The functional variation of the coincidence prob- 
ability generated by Sip will be, by definition: APi,i[ip] = 

Pi,i[4>i4> + Sip] — Philip, tp] = Pi t i[t(),ip + Sip], where 
P\ t i[ip,ip] = for identically prepared photons (coales- 
cence), as it trivially follows from Eq. (JTJ) and normal- 
ization (ip\ip) = 1. A straightforward calculation from 
Eq. ((J) yields 



APi,iM = tt 



1- a 



2 1 + A (a + a* 



(2) 



where A 2 = (Sip, Sip), a = (ip,Sip)/A and (tp,ip) = 1. 
This result is exact and perfectly general, we have not 
made hypotheses on either the magnitude or the kind of 
the functional variation Sip. In fact, we have solely used 
the basic properties of the scalar product in a complex 
linear space 5£ . 

Equation ([2]) may be further developed by assuming 
that the functional deviation Sip is generated by the vari- 
ation of a single DOF, represented by the real parameter 
/, in such a way that ip A = ip(f) and ip B = ip(f + Sf), 
with | Sf | <C l/l- For example, if / = v, the photon 
at input port A has central frequency v A = v and the 
one entering port B has central frequency v B = v + Sv. 
Defining Sip = ip(f + Sf) — ip(f) permits to express ip B 

as above: ip B = (ip + Sip) /\\ip + Sip\\ = ip + Sip. Formally 
expanding ip(f + Sf) in powers of Sf as 



iP(f + Sf) ~ iP{f) + il/ {f) Sf + ip"(f) Sf 2 /2 



(3) 



with ip'(f) = 8iP{f)/df, ij/'if) = d 2 iP{f)/df 2 , et cetera, 
we can straightforwardly obtain 



A 2 = sr 



(iP',iP')+Re(iP"^')Sf + 0(Sf 2 ), (4) 



respect to the other. The introduction of Rf[ip] has 
some remarkable consequences. Given a pair of pho- 
tons prepared in the same state ip, one can assert that 
they are maximally indistinguishable with respect to / if 
Rf[ip] < Rf[ip], V/' 7^ /. In a complementary manner, 
provided two distinct pairs of photons, the first two pho- 
tons prepared in the state ip and the second ones in the 
state (p, one can say that the photons in the first pair are 
maximally indistinguishable with respect to ip for a fixed 
/, if Rf[ip] < Rf[4>], V (p ^ ip. In summary: "all photons 
are equal but some photons are more equal than others" 
2l| . and Rf[ip] quantifies the degree of equality. 

This result conclude the first part of this Letter. In 
the remainder, we will apply Eqs. (j3]|4]) to the realistic 
case of optical Gaussian wave packets with well-defined 
spectral, spatiotcmporal and polarization DOFs. 

Gaussian wave packets. — Consider a single-photon 
wave packet of the form: 



2 

l^)=E J d 3 kMk)al(k)\0), 



(6) 



where |0) denotes the ground state of the continuous 
Fock space and a s (k) is the operator that annihilates 
one photon from the plane wave mode e s (k) cxp (ik ■ r), 
with [a s (k),al,(k')] = S^(k - k')S ss - [Hj]. Here 
{ei(k),e 2 (k),k/\k\} denotes a right-handed orthonor- 
mal basis set attached to the wave vector k. The nor- 
malization of the state is ensured by requiring 



2 

<V#)=£ Jd 3 k\Mk)f 



= 1. 



(7) 



The spectral amplitude ip s (k) determines the shape and 
the polarization of the beam. It may be deduced by im- 
posing the quantum-classical correspondence 



\E(+\r,t)\ip) 



,(+) 



(r,t), 



(8) 



and a 

Re(z) = (z + z*)/2, VzeC. Since from Eqs. fl5J} and 
flU) it follows that AP M [V>] oc A 2 = 0(Sf 2 ), we define 
the rate of distinguishability R f [ip] of the photons with 
respect to the degree of freedom / via the relation 



[(ip,ip') + (ip,ip")Sf/2 + 0(Sf 2 )] /A, where where E^\r, t) is the positive-frequency part of the clas- 



sical field wave packet whose energy equals the mean en- 
ergy of the photon in the state \ip), and 



E {+ \r,t) 



E 



RfbP] 



d*P hl [ip(f),iP(f + Sf)] 



dp 

(ip',ip') - (ip,ip' 



d 3 k 



Sf=0 



(2^)3/2 e V 2 

x exp [i (k ■ r — tot)] , 



a s (k)e s (k) 



(9) 



(5) 



where < Rf[ip] < H^'ll) an d (ip,ip') < because 
= d(tp,ip)/df = (ip,ip') + (ip',ip). The dimension- 
less parameter Rf[ip]Sf has a straightforward physical 
meaning: it tell us how rapidly two identically prepared 
photons become distinguishable when we slightly vary, 
from ip(f) to ip(f + Sf), the state of one photon with 



with uj = c\k\ = ck, c being the speed of light in vac- 
uum and eo the vacuum permittivity The expression for 
E^~\r,t) is given by the right side of Eq. © with the 
quantum operator a s (k) replaced by the classical ampli- 
tude a s (k). Then, by substituting from Eqs. ((6]) and 
(0 into Eq. ©, one obtains ip s {k) = a s (k). The to- 
tal energy contained in such wave packet is given by 
£ = /d 3 fc^(|a 1 (fc)| 2 + |5 2 (fe)| 2 ). 
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Without loss of generality, we assume that a s (k) = 
e s (k)E(k), where E(k) and e s (k) are the scalar and 
the vector spectral amplitudes of the field, respectively. 
E(k) determines the spatial characteristics of the field, 
and £ s (fc) the polarization. Here, we consider a colli- 
mated, quasi- monochromatic wave packet, with central 
wave vector fco and central frequency uj$ — c|fco| = cfco- 
We choose a normalized Gaussian spectral amplitude 
E(k) = 7(fc — fco), where 



7(9) 



(detV) 1 / 4 



r 3/4 



• exp 



(10) 



with V^ 1 — diag (of , a 2 , of) . This choice for V yields a 
factorizable spectral amplitude 7(g) = g{qi) 3(92)3(93) 
with g(q n ) = exp[-iq n r 0n - q 2 J{2o- 2 )]/ (^^^/dn). 
Clearly, it is possible to consider a more general posi- 
tive definite symmetric matrix V that couples different 
wave vector coordinates. We will see later that such a 
coupling may have dramatic consequences upon the rate 
of distinguishability of photons. In Eq. (fT0|) the real 
vector ro — {^011^02,^03} gives the position, at time 
t = 0, of the center of the wave packet. We fix e s (k) 
assuming that the wave packet has passed across a po- 
larizer that selects a uniform field polarization parallel 
to p € C 3 and perpendicular to fco, with |p| 2 = 1 and 
(fco,p) = 0. In this case, it becomes natural to de- 
fine e,<,(fc) as the normalized projection of p upon e s (fc), 
namely 0, e,(k) = (e,(fc),p)/ y/l - |(p, fc)| 2 /fc 2 , 
with |ei(fc)| 2 + |e 2 (fc)| 2 = 1 by definition. 

The Gaussian distribution 7(fc — fco) implies that the 
wave packet is concentrated in a region of the fc-space 
of "volume" o~\o~2<Jz centered at fco- Then, the assump- 
tions of collimation and quasi-monochromaticity entail 
the constraints <r, <C fco, (* = 1,2,3). In this case, the 
total energy of the wave packet can be written as 8 — 
J d 3 khoj\j{k-k )\ 2 ~ hiu Ql where / d 3 fc |7(fc - fc )| 2 = 1 
by definition. 

For quasi-monochromatic and collimated beams the 
spectral amplitude ip a (k) — e s (k)j(k — fc ) contains 
(3 + 3) + 3 + 3 = 12 independent real parameters cor- 
responding to the (spectral) © spatial © polarization 
DOFs: (fc ©{0-1,0-2, 0-3}) © r © {p e C 3 : \p\ 2 = 
1 A (fco, p) — 0}. Note that the central frequency wo is not 
an additional independent parameter, since ujq — c\ka\. 
Thus, as we claimed in the Introduction, a single photon 
is not a so "elementary" particle as it is usually believed. 
In fact, we have at our disposal as many 12 (actually 15 
if we consider a non-diagonal symmetric V), independent 
parameters to characterize the quantum state of the pho- 
ton. Each of these parameters can be taken as the vari- 
able / to evaluate the rate of distinguishability Rf[ip]. 
This calculation will be the goal of the remainder. 

Rate of distinguishability. — Using rather standard 
methods of calculation [2^, [26[ , it is not difficult to show 
that the coincidence probability (JT|) can be expressed in 



terms of the spectral amplitudes 4>f{k) and ipf(k) of the 
input photons as 



1 - 



]T /d 3 fc^(fc)^f*(fc) 



(ii) 



where fc has components {—^1,^2,^3}- This change of 
sign in the 1-coordinate is due to the parity inversion 
occurring by reflection at the BS. Hereafter, we assume 
ijjf(k) = ip s (kj) and V-f (fc) = i>s(k,f+6f). Moreover, 
for concreteness, we choose the 3-axis of the Cartesian 
reference frame directed along fco, namely fco = {0, 0, fco}. 

The explicit values of i?/, calculated from Eq. ([5]), are 
given in Table I below, for spectral and spatial DOFs: 
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TABLE I. Rate of distinguishability for several spectral and 
spatial degrees of freedom of the photons, with n = 1, 2, 3. 

A remarkable consequence from Table I, is that for the 
complementary position/ wave- vector variables, the fol- 



27] 



lowing Fourier-transform equality holds 

fifeo„^r „ = l/2, Vn = 1,2,3 



(12) 



Table I furnishes some valuable information. Consider, 
for example, the last column: it shows that i? ro „^ r 0n is 
equal to the ratio between the variation Sr 0n and the 1 /e 
width V2/o-„ of the photon wave packet in configura- 
tion space. This is in agreement with intuition: imagine 
the cross-section of each photon wave packet as a disc 
of radius \/2/o~ n . Starting from an initial condition of 
perfect superposition between the two discs, suppose to 
shift one disc with respect to the other by the amount 
Sron- Now, if Sron *C \/2/o n the two discs have still a 
large superposition and the two photons remain largely 
indistinguishable. Vice versa, if 5ro n ~ \f2/a n the two 
discs separate completely and the superposition drops to 
zero. In this case the photons become "quickly" distin- 
guishable. Analogous reasonings may be reproduced for 
the others DOFs. 

Next, we consider the case of a non-factorable spectral 
amplitude, which couples wave vector coordinates 1 and 
2. Equation (flU|) still holds, but now V has elements 
V nm = Snm/a 2 , - (S n iS m2 + S n2 6 m i)/a nm , where the real 
parameter 012 establishes the coupling, with a 2 2 > o~\o~\ 
as required by positive definiteness of V. A straightfor- 
ward calculation furnishes 



V2R a „ = 



1 1 



Or, 1 



(n = l,2), 



(13) 



where p = 01O2/0-12, with \p\ < 1. If p = (un- 
coupled DOFs) we recover the results of Table I. Vice 
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versa, for increasing coupling one has p 1 and 

R an grows unboundedly. This result is of particular rel- 
evance to experimentalists: it tells us that wave packets 
whose cross section has the shape of an ellipse whose 
either major or minor axis does not lay on the plane 
of incidence (see Fig. 1), are much more sensitive to 
mode-mismatch than cylindrically symmetric wave pack- 
ets. This fact strongly degrades photon indistinguisha- 
bility and should be avoided, for example, in coalescence 
experiments [jl. [To|. 

Finally, we examine the polarization DOFs of the two 
photons. Let us parameterize p A and p B as p x = 
{cos?9 A exp (i<Pi) , s'md x exp (i<p x ) > 0}, with A = A, B 
and \p A \ = \p B \ = 1. The results, as expansions in pow- 
ers of <j\ , <T2 , are 



of distinguishability with respect to /, /', 



2/cq 



■ cos(2tf) 



Ri 



sin 2 (2tf) 



2kg 



■ cos(2i?) 



(14a) 
(14b) 



with n — 1,2. Unlike in the spectral and scalar cases, 
R$ and R Vn are dimensionless quantities. From a phys- 
ical point of view, this means that there is not a nat- 
ural scale for the variation of the polarization DOFs. 
From Eqs. fj 14[) we see that for astigmatic wave pack- 
ets, namely for o\ ^ <72, there is a coupling between 
spectral and polarization DOFs that affects both R$ and 
Rtp n . The term (af — <j\)/{2kg) may be interpreted as 
a manifestation of the unavoidable spin-orbit coupling 
occurring in transverse electromagnetic fields [2^]. In 
addition, its absolute value furnishes the visibility of 
the coincidence fringes [29]. Equation (|14bp shows that 
R Vn cx sin(2t?). As a consequence, for linearly polarized 
states with 2d = 0, ±tt, ±27r, • • • , the phase is not a rele- 
vant DOF and R Vn = 0. 

Conclusions. — In this Letter we have studied the is- 
sue of particle indistinguishability from an operational 
point of view. In this way, we could circumvent most of 
the pitfalls usually encountered in this subtle topic [lrQ . 
In practice, from the analysis of a standard two-photon 
interference experiment we could introduce a useful pa- 
rameter, the rate of distinguishability, which furnishes a 
quantitative measure of the indistinguishability of pho- 
tons prepared in a given state. Being operationally de- 
fined, this parameter is susceptible of experimental in- 
vestigation via the measurement of the two-photon coin- 
cidence probability. 

Our main results are summarized by Eqs. (|2l5ll3p . 
and Table I. In particular, Eq. (HU1) quantifies the degra- 
dation of photon distinguishability due to coupling be- 
tween different DOFs. As a final remark, we stress that 
we have defined the rate of distinguishability considering 
the variation of a single DOF of the photons. However, 
one could also consider simultaneous variations of two or 
more DOFs and introduce the corresponding cross-rates 
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